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ABSTRACT 

We have developed a method for recovering polarization structures from the 
NRAO Very Large Array Sky Survey (NVSS) on larger angular scales than the 
nominal 15 arc minute survey limit. The technique depends on the existence 
of smaller scale fluctuations in polarization angle, to which the interferometer is 
sensitive, while the undetected total intensity of the structures can be arbitrar- 
ily large. We recover the large scale structure of the polarized Milky Way, as 
seen in single dish surveys, as well as a wide variety of smaller scale galactic and 
extragalactic features. We present a brief discussion of the uncertainties and lim- 
itations of the reprocessed NVSS polarization survey, a comparison of single-dish 
and NVSS results, and a sampling of the new polarization structures. We show 
a companion feature 1.8 Mpc outside of Abell cluster 3744, apparent Mpc-scale 
extensions to the tailed radio galaxy 3C31, a possible new giant galactic loop, and 
a new bright polarized patch in supernova remnant CTAl. We note that there is 
little quantitative information from these detections, and followup investigations 
would be necessary to measure reliable polarized fluxes and position angles. Some 
of the new features discovered in this NVSS reanalysis could provide a foreground 
for CMB polarization studies, but the internal foreground modeling for the next 
generation of experiments should have no difficulty accounting for them. 

Subject headings: galaxies: active - galaxies: clusters: general - Galaxy: structure 
- intergalactic medium - ISM: magnetic fields - supernova remnants - techniques: 
polarimetric 



Introduction 



Polarized synchrotron radiation is a powerful diagnostic tool for astrophysical relativistic 
plasmas. It provides information on the field direction, degree of order, and in some special 
situations, the actual field strength. Polarization's diagnostic power (and its sometimes 
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ambiguous interpretation) for optically thin sources arise from the combination of two factors 
- its vector nature, allowing for constructive/destructive interference, and the effects of 
radiative transfer, primarily Faraday rotation, along the line of sight. 

There is a long history to the use of radio polarizations for structural studies. In our 
Galaxy, polarized emission allowed the mapping of the la rge scale structure of the magnetic 
field including very large angle structures such as Loop 1 (jBerkhuijsen II1971I ). More recently, 
higher resolution images have u ncovered structures on a wide range of scales, including 



the in fluence of Faraday screens (ITaylor et al. 1120031 ; IWieringa et al. 1 1 19931 : iGaensler et al. 



200ll ). In addition, discrete objects such as supernova remnants and super sh ells are highly 



polarized, allowing study of their magnetic field structu res and dynamics (iKundu 1 11970 
Milne fc Dickel 1975 ; Jones et al. 12003 ; West et al. 12007 ). Many extragalactic radio sources 



were f ound to be significantly polarized, fro m the pc-scale relativistic jets of blazars and other 



AGN flWard 



1971 



Rudnick et al. Ill985l ) to the diffuse lobes of radio galaxies on scales of 



100s of kpc (jStull et al. Ill975l ). In addition to its propitious qualities, polar ized synchrotron 



emiss ion is also important as a foreground for CMB polarization studies (ITegmark et al. 
200(]h . 



Polarized objects like th e Mpc-scale peripheral relics or gischt (e.g.. lHanisch. Strom fc Jaffe 



1985 



Kempner et al. 112004) around clusters of galaxies, presumably due to accretion shocks 



Ensslin et al. 19981 ) were the initial motivation for the current study The current work was 



motivated by a desire to use polarization as a probe to push below the limits of confusion 
( IRudnick 112004 ) to search for signatures of infall into and alo ng the large scale filaments that 
are part of cosmic structure formation (jMiniati et al. Il200ll ). 



Although the Stokes parameters Q, U, and I for linearly polarized emission contain a 
great deal of information, the way in which they are processed and communicated necessar- 
ily destroys some of that information. In particular, whether or not a vector background 
subtraction (in Q and U) or spatial smoothing is done before calculation of the polarized 
intensity P = [Q 2 + £/ 2 ] - 5 can cause features to appear or disappear in an image. Thus, it 
is often possible to re-process Stokes images and derive further information than emerged 
from their initial analysis. It was this type of thinking that led us to reconsider the pol ariza- 



tion information present in the NRAO VLA Sky Survey (NVSS, ICondon et al. Ill998l ). We 



note that the proper way to recover diffuse polarization information unconfused by smaller 
sources is to combine single dish and interferometer data. However, since there are angular 
scales and rotation measure ranges now reachable by the "allsky" NVSS that are unlikely 
to be duplicated in the near future, it is worth examining what information can be derived 
from this survey. 



The NVSS is a 1.4 GHz survey (combining data from 1364.9 and 1435.1 MHz) conducted 
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from 1993-1997 in the D and DnC configurations (lowest angular resolution) of the VLA, 
covering the 82% of the sky which is visible from its latitude. It is reproduced as 2326 4° x 4° 
images in Stokes parameters I, Q, and U at a resolution of 45". The rms brightness is ~0.45 
mJy/beam (0.14 K), in I, with typical rms values of ~0.29 mJy/beam (0.09 K) in Q and 
U at high galactic latitudes. Faraday rotation between the two NVSS bands eliminates any 
polarized emission in the images with |RM| > 340 rad/m 2 . 

The information on large scale structures in the NVSS is limited by two effects. First, 
the snapshot observations and smallest projected spacings of ~37 m produce grating lobes 
at ~18', effectively eliminating the sensitivity to sources larger than that scale. In addition, 
the short spacings were further weighted down in the mapping procedure, to remove the 
large scale "pedestal" in the synthesized beam. Although no detailed study has been done, 
the NVSS is generally recognized to have little sensitivity to structures >15'. For most 
extragalactic studies, this provides a great benefit both in total intensity and polarization 
because it removes the effects of the strong large-scale galactic emission. However, there 
are many interesting galactic and extragalactic structures on intermediate scales that do not 
appear in the original survey. 

It is possible to recover some of this intermediate scale structure from the NVSS linear 
polarization images. If a source has structure in Q and U on angular scales to which the 
interferometer is sensitive, then it will be detected even if the underlying total intensity (I) 
is much smoother and undetectable by the interferometer. Q and U often have smaller-scale 
fluctuations than total intensity because they are sensitive to field disorder, magnetic field 
direction, and Faraday rotation, none of which affects the total intensity. The recovery of 
such polarized features is the purpose of our reanalysis; the NVSS "all sky" images in total 
intensity and polarization at 800" resolution are shown in Figure [lj dramatically illustrating 
the larger scale information present in the latter. The processing scheme for these images is 
discussed in detail below. 



1.1. Polarization Image Preview 



We briefly look at what types of features and artifacts are visible in Figure [T] and the 
additional closer look in Figure [2] (center image) at a ~100° region centered around the 
Perseus arm. First, we see that the polariz a tion i mag es trace the la r ge sca le emission of 



the Galaxy, as seen, e.g., in iReich Sz Reich I Jl986h or lHaslam et al. I (119991 ). Superposed 



on the galactic structure are a network of stripes following lines of constant right ascension. 
These lines represent series of scans from the NVSS survey where the residual instrumental 
polarization was slightly higher than average. There are also lines of black 'dots', each of 
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them one approximately one primary beam (34') across, where the residual instrumental 
polarizations were high enough to require flagging in the NVSS survey (W. Cotton, private 
communication, 2008). At the highest declinations, the residual instrumental polarization is 
high at all right ascensions, creating a series of rings around the north celestial pole. 

A large number of small sources (at 800" resolution) are also seen in these images. These 
represent both real polarized sources and those caused by residual instrumental polarizations 
from the very brightest objects in total intensity. To illustrate the instrumental issues from 
strong point sources, we consider the case of 3C84; it has an NVSS flux of 22.15 Jy, and an ob- 
served polarized flux of 1.9 mJy (0.08%), showing that the residual instrumental polarization 
is quite low at full resolution. However, larger sources such as Cassiopeia A (~5' diameter) 
seen in the center of Figure [2] are incompletely sampled in the NVSS, producing sidelobe 
structure that puts power into the polarization images. Out of Cas A's total flux of ~2500 Jy, 
the total polarized flux in Figure [2] is ~45 Jy, or ~2%, sprea d over Although there may 



be so me small real integrated polarization at 1.4 GHz (see I Anderson. Keohane fc Rudnick 



19951 ). the 1° extent shows that the bulk of the contribution here is instrumental. Thus, very 
strong, extended sources cannot be reliably studied through this NVSS reprocessing. We 
consider this effect further below in our discussion of emission around 3C31. 

By contrast with Cas A, there appears to be little or no effect due to the strong total 
intensity emission from the Galaxy , e.g., on scales of degrees or larger. This potentially could 
have been a problem, as lStil et al. I (120061 ) show that the noise in VLA visibility data is related 



to the value of T sys , which includes contributions both from the sky brightness and spillover 
radiation from the ground. However, even in the brightest regions of the galactic plane, we 
do not appear to have preferentially high polarization; the brightest region in Figure [2j for 
example, is in the direction of the Cygnus Arm, and is indicated by a black circle. Withi n that 
region, the peak (mean) brightness is ~35 K (~7 K) from the iReich &: Reich I (119861 ) map, 
while we observe a mean polarized brightness of ~20 mK above the background. Instead of 
being higher in this bright region, the NVSS polarized flux is strongly anti- correlated with 
the total intensity e mission, likely due to d epolarization fr o m the extensive ionized gas seen 
over this same area (IHaffner et al. II2003I ). IStil fc Taylor I (120071 ) also find a major drop in 
the number of polarized NVSS sources in this area, likely due to depolarization between the 
two NVSS bands. 



A subtle problem with NVSS polarizations has been described by lBatttye. Browne fc Jackson 



(120081 ). largely resulting from clean biases and small (//Jy) offsets in Q and U. Both of these 
create problems far below the noise in the NVSS polarization images, and even further below 
the residual instrumental polarization artifacts that are prominent in Figures [1] and [2], and 
do not affect the current work. 
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2. Map Analysis 

The key element of our processing is the convolution to larger angular scales of the 
polarized intensity maps, instead of the typical procedure of pre-convolving Q and U, and 
then calculating P. Such pre- convolution is generally preferred due to its better noise prop- 
erties, as well as its preservation of the polarization angle (x = 0.5 x tan^ 1 {Q/U) ) on the 
convolved scales. Pre-convolution is also equivalent to observations made with a large syn- 
thesized beam or single-dish observation. However, since there is a maximum angular scale 
to which interferometers are sensitive, pre-convolution at that maximum scale or beyond 
would contain no real signal, only noise. 

Once P is calculated from Q and U images at some resolution, further smoothing (which 
we term post- convolution here) can be performed and arbitrarily large structures can be 
detected. Similarly, if a large total intensity structure were not smooth, but made up of a 
collection of small-scale features, it would also be detectable beyond the nominal maximum 
angular scale. This process is diagrammed in Figure [3J It shows the production of two 
maps, P 8 oo an d Psoo/ > constructed as described above by pre-convolution followed by post- 
convolution at a resolution of 800" . 

The rms scatter in P 800 and Psoo/ varies widely across the sky, and contains contribu- 
tions from three factors - the random noise from the receivers, the actual polarized signal 
in each region, and the residual instrumental polarization which is most easily evident in 
the strong declination stripes visible in Figure HJ During the NVSS observations, the focus 
was on the reliability of the total intensities, and changes in the instrumental polarization 
were not closely monitored (B. Cotton, private communication). At times when certain 
declination strips were being observed, larger than normal instrumental effects would arise, 
e.g., from phase jumps between the right and left hand receivers on the reference antenna 
used for polarization calibration. We illustrate more quantitatively in Figure [8] how the rms 
fluctuations in P 800 vary across the sky at high galactic latitudes. We identify the high rms 
scatter region associated with the actual signal related to the North Polar Spur, and an adja- 
cent region where the residual instrumental polarization contributions were small. These can 
also be seen visually in Figure HJ In the pr esence of random noise alone, we could have cor- 



recte d for the bias in polarized intensity (IWardle Kronberg I Il974l ; ISimmons fc Stewart 



19851 ). Given the dominance of variable residual instrumental polarization contributions, 
however, it is not possible to conduct an automated processing scheme for separating the 
signal/noise/instrumental contributions to the power. Therefore, we look at the noise charac- 
teristics for each individual source of interest, and describe the observed value in P 80 o above 
the local background and rms to estimate its significance. 



To remove some instrumental artifacts (and a first order subtraction of the polarization 
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bias) we will also be making use of a high-pass filtered version of Psoo , which we label Psoo/ 
. The high-pass filter consisted of subtracting from each pixel the median of Pgoo over a box 
4° in declination (the size of an individual NVSS image field) by 6' in right ascension. 

We estimated the sensitivity to diffuse polarized structures after pre- and post-convolution 
using a series of simulated polarized signals. The initial construction of the simulated images 
(Signal+ Noise) is shown in Figure HI the processing of these images to measure sensitivities 
is shown in Figure One of the processed (Signal+Noise) images is shown in Figure El In 
order to capture all the spatial characteristics of the NVSS images, we used actual images 
from high galactic latitude fields for both the simulated "signal" and "noise" . We constructed 
the background "noise" image by first taking a single pair of 8.5°x 8.5° Q and U images 
from the survey, and at all locations where |Q| and |U| were greater than 2 mJy/beam, 
replaced these values with zero. This effectively eliminated almost all signals from compact 
sources in the Q and U images. The remaining large-scale real and instrumental fluctuations 
in these images thus simulate the actual background against which we are trying to detect 
signals. The rms "noise" distribution, which includes point source, galactic and instrumen- 
tal fluctuations, is shown in Figured for all 4° fields in the NVSS. It can be very roughly 
described as a Rayleigh distribution with an rms value of ~6.5 mJy/(800" beam), slightly 
larger than the expected system noise value of 5.2 mJy/ (800" beam) = 0.29 mJy/(45" beam) 
x([800"/45"] a5 ). At the full resolution of the NVSS, the systematic real and instrumental 
effects contribute little power to the Q and U maps. The systematic contributions become 
much more important as the images are convolved to lower resolution. 

To create our simulated signal, we took a separate NVSS 1° field and convolved the Q 
and U images with an 8' Gaussian, multiplied these images by a constant and added them 
respectively to the background noise Q and U images. The simulated signal thus consists 
of a 1° patch with structures in Q and U on the scale of 8', added to the 8.5° background 
(noise) field at the full NVSS resolution of 45". 

We then pre-convolved these images to various scales, converted to P (first option in 
Figure [5]) and measured the means in the central 1° region (P S i g , including the local noise) 
and outside the central 1° (J 3 noise)- We also measured the rms fluctuations, P Tms in the noise 
regions. Finally, we calculate 

Signal : Noise = [P sig - P noise ]/ 'P rms - 

The resulting signabnoise as a function of pre- convolution size is shown in Figure [7J The 
detectability rises dramatically with convolution, as expected, since the "signal" was forced 
to have a scale of 8'. At larger pre-convolution values, the signabnoise falls as the signal 
itself starts to be vector-averaged away and diluted. 
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We also show the behavior of the signal:noise using post-convolution. For these experi- 
ments, we added the "signal" image (with its 8' scale) to the full resolution Q and U "noise" 
images, formed P, and then post-convolved the result to various scales (second option in 
Figure E]). Looking again at the signal:noise for various convolutions, we see that the de- 
tectability of the signal is less than in the vector-averaged (pre-convolution) case because the 
signals are not being averaged coherently. At the larger convolution sizes, the signal:noise 
drops even further, mostly due to the increased noise power on large scales in the NVSS, from 
both galactic polarization structure and residual instrumental polarization. We note that 
the above signal detectability experiment is a single case using an arbitrarily constructed 
source. The detectability of sources with specific polarization structures would have to be 
investigated on an individual basis. 

It is also possible adopt a hybrid approach and pre-convolve an image to some scale 
(Opre), and then post-convolve it to larger scales (6 pos t)- The advantage to this approach 
is that if an astronomical source has structure in Q and U on scales ~ source , then the 
signahnoise will be enhanced by using 9 pre ~ 6 SOUTCe . This pre-convolution has a second 
advantage for our survey - it enhances our search for new larger-scale features by enhancing 
them relative to smaller sources that are more likely to be already known. We therefore 
adopt, for this initial presentation 9 pre = 240" and 9 post = 800", and call the polarized 
intensity Psoo/ • It is important to note that while our Psoo/ maps with 6 post = 800" have 
a resolution very similar to those of single-dish surveys, sources with 6 SOUTCe ~ 240" will be 
detected by us, but not with the single dish. 

The results presented here are almost entirely qualitative, for a number of reasons. First, 
they measure only the polarized "power" in Q and U from the poorly sampled low spatial 
frequencies in the survey up to the pre-convolution scale size. Thus, a source one degree 
across with uniform Q,U would be invisible, while another source with variations on, e.g., 
4 arcminute scales would appear strongly, even if they had the same mean sqrt(Q 2 + U 2 ). 
The spatially dependent power in Q and U may also depend on the foreground rotation 
measure, independent of the intrinsic strength of the source. Quantitative information could 
be derived only by constructing a detailed polarization model of a source and its foreground 
rotation measure structure, and then propogating that model through simulated observations 
and processing. We do not attempt this here. 

3. Results 

We begin with a discussion of two large fields in the galactic plane, to illustrate the 
differences between structures seen in total intensity and/or single-dish polarization maps. 
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We also look at several new diffuse polarized structures, galactic and extragalactic, found 
through our new processing. A paper identifying a large number of new and often unidentified 
sources is in preparation. 



3.1. Cygnus-Perseus region 



This region is shown in Figure [2] and we have previously discussed some of its instru- 
mental and other features. Here, we point out the filamentary feature extending north for 
about 40° from the galactic plane in the Cygnus region (1/7=90°) with a less-well defined 
counterpa rt to the south. The fil ament is also present in the total intensity Bonn image at 36' 
resolution iReich fc Reich I (119861 ). especially when filtered to remov e the smooth background 
(using the multi-resolution filtering technique of iRudnick I (120021 ) with a box size of 7° in 
long itude by 1° in lat itude). This feature is usually identified as the western portion of Loop 
III (ISpoelstra 1119721 ), but another pos s ibility is suggested in Figure [HI The 36' resolution 
Stokes Q image from IWolleben et al. I (120061 ) is shown here, with the suggestion of a new 
loop that extends from the Cygnus region, curves up and over the galactic pole, and returns 
back towards the galactic plane approximately 180° away. With a nominal center at lj/=0°, 
b//=20°, and a radius of 75 ° (M. Wolleben, p rivate communication), it is approximately con- 
centric with Loop I, which Wolleben (j2007l ) has recently modeled in terms of two ~100 pc 



radius synchrotron emitting shells in which we are embedded. This new suggested loop needs 
further study through other ISM tracers, and could potentially revise our understanding of 
our local environment. 



3.2. Galactic Center region 

A 50° field near the galactic center (Figure [TU1) further illustrates the similarities and 
differences between single dish polarization measurements and the NVSS polarization re- 
processing at the same wavelength and angular resolution. We pre-convolved the NVSS Q, 
U images with 240" and post-convolved the P images with 36'. In the NVSS image, the 
vertical bands are artifacts from small differences in the residual instrumental polarization 
at different times, as observations were made along lines of constant right ascension. This is 
also likely responsible for some of the fine-scale mottling apparent along the vertical bands, 
but it is not possible for us to separate these from fine scale polarization structures. We 
believe that the rest of the structures visible in this NVSS image are "real" in the sense that 
they reflect the actual signals from the sky. Their interpretation, however, is quite different 
from those of single dish images, as described below. 
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As is well known, polarization and total intensity structures are very different, even 
when th ey are both from single -dish measurements. The image from the DRAO polarization 
survey JWolleben et al. I bood ) shows a bright patch in the upper right, with a smoothly 



varying position angle (not shown) over scales of degrees. The NVSS is not sensitive to Q 
and U variations on such large angular scales, and the feature is not visible in our images. 
High rotation measures (>340 rad/m 2 ) would also cause a degradation of the NVSS signal 
compared to the DRAO image. Along the galactic plane, running from upper left to lower 
right, the DRAO image shows an unresolved narrow bright band bordered by two dark 
stripes. The dark stripes are beam-depolarized due to a rapid switch in polarization angle 
across these lines. In the NVSS data, these same patterns in the sky create a different 
response; a broken thin strip of polarized emission along the plane is seen from the regions 
where the angle is changing rapidly, creating a detectable interferometer signal in Q and U. 
However, farther from the plane, the relatively constant polarization angles from the DRAO 
position angle image (not shown here) lead to a broad dark region in the NVSS image. The 
numerous depolarization filaments seen elsewhere in the DRAO image are again not seen in 
the NVSS; rapid changes in polarization angle provide an NVSS signal at high resolution, 
and there is no depolarization introd uced by the po s t-con volution. An alternative probe 



of depolarizing regions is presented by IStil fc Taylor I (120071 ). who study the polarization of 



compact sources in the NVSS. They find regions ~10° across that cause a reduction in the 
number of polarized NVSS sources, due to both intervening H II regions and diffuse galactic 
structures. 

In Figure [TU], we also point out the bright emission from 3C353 in the NVSS image. This 
radio galaxy has very strong polarized (Q and U) emission when convolved to the 240" scale, 
which is then carried forward through the post-convolution to 36'. However, the averaging 
of Q and U themselves over 36' in the DRAO image makes 3C353 undetectable against the 
galactic background. 

The very luminous HII region M17 is bright in both the DRAO and NVSS images. The 
NVSS apparent extended polarized flux comes from strong sidelobes, as can be seen in the 
full resolution images; it is likely that the DRAO polarization is also due to instrumental 
polarization. W44 is detected strongly in DRAO and weakly in NVSS. The full resolution 
NVSS polarization image (not shown here) shows a structure which is largely unrelated to 
the total intensity structure of W44, so the physical origin of the polarization seen at 36' 
resolution is unclear. 
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3.3. The Extragalactic Sky 



Away from the galactic plane, we get a better view of more compact polarization features 
in the NVSS, as seen in Figure [TT1 This is an ~ 44° x 17° strip of the sky centered around 
RA, Dec 11.8 h, 34.2°, with Psoo/ i n re d. The green image shows the ROSAT All-Sky Survey 
(IRASS 1 120001 ) convolved to 800". The blue image is the total intensity NVSS image, also 
convolved to 800". Note that at this high galactic latitude (35° < bn < 90°), and with the 
additional filtering, very little galactic structure is visible in Psoo/ • A large number of bright 
compact features are visible, appearing in red (blue) for high (low) fractional polarization. 
Most of these are also seen in the full resolution NVSS images. 

A number of more extended polarization features are also visible in Figure [TT], many 
of which are well-known sources. In the southeast, the bright green structure is the X- 
ray emission from the Coma cluster of galaxies, w ith an extension to the southwest from 
an inf ailing sub cluster (IFeretti fc Newmann II2006I ). At the southwest X-ray terminus is a 
transverse polarized radio structure with brightness ~20 mJy/80 0" beam (20 mK ) abov e 



the background. This is the well-studied "relic" first detec ted by IJaffe fc Rudnick 



and then studied in total intensity and po larization, e.g., by lHanisch. Strom fc Jaffe 



1979), 



1985|); 



Giovannini. Feretti fc Stanghellini I (11991I ). Note the absence of blue, total intensity, emission 
from this polarized structure; with an angular extent of ~1°, the total intensity is not 
detected in the NVSS survey. Similarly, the diffuse polarized emission from the giant radio 
galaxies 3C236 and B2 1321+31 is easily visible, while only their more compact features are 
seen in total intensity in the NVSS. 

Among the bright extended polarization structures are several that were previously un- 
known. Feature A is one of the brightest structures in the high latitu de sky (25 - 35 mK); 



follow -up observations with the Westerbork Synthesis Radio Telescope ( iBrown and Rudnick 



20081 ) suggest that it is a galactic Faraday system. Such structures have no intrinsic syn- 
chrotron emissivity, but appear bright to the interfero meter because they produc e small-scale 
variations in Q and U from the galactic background (IWolleben fc Reich 1 120041 ) . Feature B, 
and others like it in the image, are examples of single NVSS pointings with high residual 
instrumental polarization; they are recognizable by their circular shape and sizes comparable 
to the primary beam (~30'). We have not yet found a reliable way to eliminate instrumental 
problems at the lowest levels seen in Figure [TTJ especially when they span several pointings. 
In the following, we therefore present examples of new extended polarization features only 
with high signal:noise, to illustrate the potential power of this NVSS reprocessing technique. 
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3.4. Abell 3744 



A network of tailed radio sources is seen in this cluster (jMarvel. Shukla fc Rhee Ill999l ). 
which might be responsible for the bright polarized emission centered on the cluster seen in 
Figure [121 At a redshi f t of .0381, Abell 3744 is listed as a member of supercluster number 
180 by lEinasto et al. I (120011 ). who group it with Abell 3733, (z=0.0382), and, apparently 
mistakenly, with the background clus ter Abell 3706 at a re dshift of ~0.1 . Abell 3744 
is a member of the REFLEX sample (IBohringer et al. 1 120041 ) with a relatively low X-ray 
luminosity of 1.6xl0 43 erg/s, integrated over 0.1 - 2.4 keV. The X-ray emission is seen only 
in the region of the cluster center. 

To the east, a polarized structure is seen 37' (1.8 Mpc) from the cluster center, with an 
extent of ~1.4 Mpc. A slice at constant declination through the peak of polarized emission 
and through the eastern structure is also shown in Figure [T2| to demonstrate the signal: noise 
of these features. The eastern structure has a peak polarized flux 40 mJy/beam above the 
background, with a background rms of 6 mJy/beam. It has a total polarized flux (Psoo/ ) 
of ~90 mJy, and no obvious total intensity counterpart. It is important to note that our 
polarized fluxes are the result of our non-standard processing procedure, and that the true 
signals measured with an interferometer plus single-dish system could be considerably higher, 
assuming no additional contributions from a nearby Faraday screen. Assuming our nominal 
fluxes and a 33% fractional polarization would yield a monoc hromatic radio lumin osity of 
10 24 W/Hz at 1.4 GHz. If this is a peripheral relic, or gischt (jKempner et al. 1120041 ). both 
its radio and associ ated X-ray cluster luminosities are significantly lower than those of most 



relic systems (e.g.. iGiovannini. Tordi &: Feretti Ill999l ). O nly one other cluste r, Abell 548b 



which is part of a very complex optical and X-ray system (IDavis et al. Ill995l ). has similarly 
low luminosities. The Abell 3744 relic would also be considerably further from the cluster 
core than is typically seen. Confirming and more detailed polarization observations of this 
system would also be most useful. 



3.5. 3C 31 



3C31 is a very well-studied wide-angle tailed radio galaxy (IFanaroff &: Rilev II1974T) with 
a tot al extent of >40' (~850 kpc at z=0.0169 with H = 70 km/s/Mpc. iKlein fc Wielebinski 
19791 ). In our polarization image, Figure [T3J, we find emission extending ~1.2° to the south 



west, where it partially merges with polarized emission from an unrelated background source. 
To the southeast, there is polarize d emission bridging to 3C34, an unrelated source in 
a z=0.69 cluster ([McCarthy I Il998l ). Formally, the statistical signifance of the polarized 
emission is high. Towards the southwest (southeast), the polarized brightness is ~15 (25) 
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mJy/800" beam above the background, with a background rms of 4 mJy/800" beam. How- 
ever, there is an additional source of contaminating polarized emission in the neighborhood 
of strong, very extended polarized sources such as 3C31. These are sidelobe structures from 
the poorly sampled short baselines, but would appear as true positive-definite signals in our 
processing. To assess the importance of this effect near 3C31, we constructed an equivalent 
'Jsoo" image by processing the total intensity image in an identical way to Psoo , substituting 
I for both Q and U. This created excess power in the vicinity of 3C31, as expected. How- 
ever, the ratio of the brightness in the southwest (southeast) polarized feature to the peak 
brightness of 3C31 itself is 0.1 (0.17) in Pgoo , whereas the brightness ratios for the same 
locations in J 8 oo are <0.01. We therefore conclude that poorly sampled sidelobes from the 
polarized emission in 3C31 is not responsible for the newly detected features. 

It is not clear whether this diffuse emission is associated with 3C31 or 3C34. The host 



galaxy of 3C31 is NGC 383, the brightest of a rich group of galaxies (IBurbidge fc Burbidge 



196ll ). many of which are seen in the area of the polarized emission to the SE of the sourc e 
(Figure [5]). These are all part of a filament of galaxies at this redshift ( jMiller et al. 1120021 ). 
which extends through the SW polarized extension over 9° to the NE to cluster Abell 26 2 



(IMoss fc Dickens 1119771 ). In the A 11.1 cm single-dish images of iKlein fc Wielebin ski 



a 10' e xtension is seen to the southwest. The low resolution 102 MHz images of lArtyukh. Ogannisyan &: Tyu 
(J1994J) also show extended emission to the southwest, with a size <50', the extent of their 
primary beam. There is possible contamination from the southwest background source. Note 
that in their Figure 1, 3C31 is the northern component of the apparent double structure; 
the southern component is 3C34. 

Assuming that the southwest extension is associated with 3C31, it would have an extent 
of ~1.5 Mpc. There is some possibility that the observed polarization structure results from 
distant sidelobes of 3C31, although processing the total intensity NVSS images in the same 
way as Q and U (i.e., forcing them to be positive and then post-convolving), did not show 
significant sidelobe contributions. Assuming the polarization structure is real, the All.l cm 
results cited above also suggest that there has been an outflow to the southwest from 3C 31. 
Lacking bright jet features, this structure would likely be a remnant of past activity, perhaps 
now being energized by large scale group processes. Indications of the dynamical state of the 
group include the offset of NGC 383 fr om the centroid of the X-ray e mission and significant 
X-ray structure towards the southwest (IKormossa fc Bohringer Ill999l ). which is detected out 
to the virial radius at ~700 kpc. 
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3.6. CTA1 



This well-known supernova remnant has a filled X-ray structure and a radio shell which 
is not visible in the NVSS total intens ity images. It has bee n mapped at 1' in polarization 
at 1.4 GHz using the DRAO telescope (IPineault et al. 1119971 ). All of the major features seen 
in the DRAO image are reproduced in our NVSS reconstruction at 800" resolution (Figure 
[T5l) . In addition, we find one diffuse bright patch at OOhlOm, 72° that is not seen in the 
DRAO map. The diffuse nature of this patch is established by the lack of smaller-scaled 
polarized features in the full resolution Q, U images (at a level of Q, U < 3 mJy/45" beam). 
Its peak brightness is 30 mJy/800" beam above the local mean, with a background rms of 
~5 mJy/800" beam. If it were completely smooth, the polarized brightness would be only 
0.2 mJy/r beam, compared to the DRAO rms value of 0.3 mJy/1' beam, so it is reasonable 
that it had not been previously detected. 



The bright patch is on the border of what IPineault et al. I (119971 ) call the "reverse shell" 
region, where the curvature of the radio structure is inverted. Their explanation for the 
shape of this region is that the shock has encountered and encircled a dense cloud. The new 
polarized patch may help define the borders of this cloud enc ounter, but the l ack of bright, 
narrow polarized or total intensity emission, as pointed out by IPineault et al. I (119971 ). is still 
a problem for this explanation. 



4. Discussion 

Polarized intensity is not a well-defined quantity when there is more than one polar- 
ized component along the same line of sight, creating both challenges and opportunities. 
Although there are heroic attempts to separate multiple components if they are separated 



in ro tation measure space (IBrentjens fc de Bruyn 1 120051 ; ISchnitzeler. Katgert fc de Bruyn 



20071 ). in general, observed polarized intensities depend on a complicated combination of 
observing frequency, resolution, and the spatial frequencies represented in the image. This 
complication also allows a variety of processing schemes to be carried out on Stokes param- 
eter images, which can result in new structures being identified. 

Our reprocessing of the NVSS survey in Stokes Q and U explores a region of observing 
frequency/resolution/spatial frequency that has not been previously studied. We therefore 
have been able to identify many new features, as well as recover structures that are invisible 
in the total intensity NVSS but seen at other telescopes. One common feature of the ex- 
tended emission features we detect is that they are of low surface brightness. Psooj images 
free of galactic emission have rms fluctuations ~2 mJy/800" beam (2 mK). A 5a detec- 
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tion would therefore correspond to ~10 mJy/800" beam or ~30 mJy/800" beam in total 
intensity, assuming a fractional polarization of |. By contrast, a single dish working at the 
same frequency and angula r resolution, would have a 5a confusion limit of ~125 mJy (e.g., 



Condon fc Broderick Ill985l ). Under favorable conditions of high fractional polarization and 



Q,U variations on scales <15', the reprocessed NVSS can then be much more sensitive to 
diffuse structures than any single-dish survey. 

All-sky surveys of diffuse polarization could thus be a powerful way to probe low density 
extragalactic regions. Assuming an optimum sensitivity of 30mJy/800" beam (total inten- 
sity), the equivalent minimum energy magnetic field for a detected synchrotron structure 
will have values ~0.2uG. T his corresponds to pressures of ~10~ 14 ' 5 in the regime of the 



WHIM (IKang et al. II2005I ). Radio observations at these low brightness levels could then 



serve to illuminate the diffuse baryon component of large scale structure, where 30-50% of 
the baryons are likely lo cated, bu t are e xcee dingly difficult to detect through their thermal 
emission. Simulations by iMiniati I (120041 ) and lPfrommer et al. I (120061 ). e. g., show that such 
radio emission is expected, driven by the magnetic field amplification and relativistic particle 
acceleration at shocks on the borders and in the interior of filaments. 

Much more sensitive polarization measurements, with significantly reduced instrumen- 
tal problems, will soon be available on the EVLA0, and a new generation of wide-field, low 
frequency instruments such as LOFAR0 and the MWA| will allow us to probe much deeper 
into the WHIM regime over large areas of the sky. Combinations of single dish and interfer- 
ometer observations may also allow such features to be seen in total intensity for individual 
regions where the large investment of observation and analysis time can be justified. 



4.1. CMB polarized foregrounds 

It is unlikely that features found through our NVSS reprocessing will affect upcoming 
CMB polarization experiments. We have identified new features on scales down to 9 pre 
= 240" . For the new synchrotron sources, forthcoming CMB polarization exper iments that 
work at these resolutions or smaller, such as Planck or EBEX JOxlev et al. 1120041 ) .H will have 



sufficient frequency coverage to model the foreground contamination from galactic or extra- 



1 http:/ /www. aoc.nrao.edu/evla 
2 http://www.lofar.org 

3 http://www. hay s t ack .mit.edu/ast/arr ays / mwa 

4 See http://lambda.gsfc.nasa.gov for a complete list of CMB polarization experiments. 
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galactic synchrotron sources internally. We also detect "pseudo-sources" due to Faraday 
screens in our galaxy. However, direct Faraday rotation of the CMB signal from these 
regions is negligible at v > 100GHz (Ad ~ 0.5 deg for a rotation measure of RM=1000 
rad/m 2 ). At the same time, the Farada y screen features could indicate the presence of 
undetected HII regions (ISun et al. 1 120071 ). The bremsstrahlung emission from HII regions 
can be polarized at the 10% level due to Thomson scattering at the edges of the clouds 



Keating et al. I (119981 ). but the polarized brightness should be at least an order of ma g nitud e 



less than that of synchrotron emission at frequencies above 10 GHz iBennett et al. I (119921 ). 



Such sig nals are typically not modeled as a p olarized foreground component in the WMAP 



analysis (IPage et al. 1120071 : iGold et al. 1120081 ) 



5. Conclusions 



A reprocessing of the NVSS polarization images has allowed the recovery of diffuse struc- 
tures on large angular scales. The details of the processing allow one to tailor the sensitivity 
to particular angular scales of interest. At present, residual instrumental polarization vari- 
ations across the sky are a key limiting factor. A variety of new galactic and extragalactic 
sources have already been identified, with a more comprehensive census underway. While a 
better recovery of diffuse polarization structures is possible, e.g., by combining single dish 
and interferometer measurements, the next generation of radio telescopes such as LOFAR, 
the MWA and EVLA will also be able to exploit the processing technique introduced here 
to provide probes, e.g., of the relativistic plasmas associated with the elusive WHIM. 
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Fig. 1. — All sky AITOFF projections in galactic coordinates of the NVSS, at a resolution 
of 800". The top image shows the total intensity, with polarized intensity on the bottom. 
The obvious striping in the polarized intensity, mostly along lines of declination, are due 
to variations in the residual instrumental polarization during the survey. Occasional small 
black regions is where the residual instrumental polarization was especially high, and the 
polarization data were flagged for that individual pointing. 
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Fig. 2. — Total and polarized intensity around the Perseus-Cygnus region of the galactic 
plane. Galactic coordinates are indicated on the left image. "X" marks the position of the 
celestial pole; circles around the pole visible in the middle im age show lines o f con stant 
declination. Left: 21 cm total intensity from Bonn 25m survey iReich fc Reich I (119861 ). 36' 
resolution. Center: Polarized intensity from NVSS at 800" resolution, with no filtering. 
Right: Filtered version of map on left, as described in text. The black circle marks the 
Cygnus arm at ljj=90°; the arrow indicates the base of what is normally called "Loop III" 
and is discussed further in the text and in Figure [9j 
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Fig. 3. — Flow chart for construction of polarization images P 80 o and Psoo/- 
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Fig. 4. — Flow chart for construction of simulated images (signal+noise) for sensitivity 
experiments. The NVSS Q and U images were first clipped to remove the effects of strong, 
compact polarized sources. 
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Fig. 5. — Flow chart for two options of processing (Signal+Noise) images for sensitivity 
calculations. 
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Fig. 6. — Example of simulated polarization data used to calculate sensitivity. The central 
1 square degree contains the simulated signal, a random pattern of polarized flux with a 
characteristic scale of 8'. In this particular example, the image was "pre-convolved" by 4' 
and "post-convolved" by 8'. 
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Fig. 7. — Left: Histograms of the rms noise in each 4° field in the NVSS. Right: signahnoise 
from simulated polarization measurements, as described in the text. The black line indicates 
"pre-convolution" of the Q and U images; the grey (lower) line indicates "post-convolution" 
of the P image. 
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Fig. 8. — RMS scatter in P 800 in strip from 53° < b n < 67° as a function of galactic longi- 
tude. Inidividual spikes and the large scatter are due to variations in residual instrumental 
polarization. A region of low instrumental effects is indicated, along with the increased rms 
power due to the North Polar Spur. 
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Fig. 9. — DRAO Stokes Q image (jWolleben et al. 1120061 ) in celestial coordinates, centered 
at RA„ Dec. ~ 15h, 30° , and approximately 150 ° in RA by 110° in Dec. "L" marks the 
base in the galactic plane of the suggested new loop, at \ H =9>2°. while the arrow shows the 
position of the galactic center, just below the edge of the image, so the galactic plane runs 
between these two. X marks the galactic pole, and NPS denotes the North Polar Spur. A 
curved line is drawn outside of the possible new loop. 




Fig. 10. — Left: Polarized intensity from DRAO survey (iWolle ben et al. 1 120061 ) . Right: 
Polarized intensity from NVSS. Both images have a resolution of 36'. The field is centered 
at l,b = 26.5°, -1.1° and is 50 degrees in width. 
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Fig. 11. — Image of approximately 44°xl7° strip in celestial coordinates, centered around 
11.8h,34.2° Red is Pgoof , green shows the broadband X-ray emission from ROSAT, convolved 
to 800", and blue is the total intensity emission from NVSS, convolved to 800". 
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Fig. 12. — Abell 3744, contours and slice of NVSS polarized intensity at resolution of 800" . 
The top image is overlaid by the NVSS total intensity at 45" resolution. Contour levels are 
at 0.01 x (4, 5, 6, 8, 10) Jy/beam. The total field is 95" in declination by 115" in right 
ascension. The slice is taken in right ascension through the peak and extends over 5° to 
show the signahnoise of the detected features. 
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Fig. 13. — 5.3° field with 3C31 (upper left) and NGC 315 (lower right). Psoo/ is in red. 
The brightness of the southwest extension of 3C31, ignoring the region of the background 
double, is ~12 mJy/800" beam. The vertical and horizontal lines indicate slight differences in 
background removal from the original 4° NVSS fields. The green image is the full resolution 
WENSS survey image at 330 MHz. 




Fig. 14. — Contours of polarized emission (Psoo/ ) in a 2° field around 3C31, at levels of (6, 
13, 28, 38, 51, 76, 100, 126 and 140) mJy/800" beam. Small cir cles indicate the pos itions of 
galaxies at redshifts of 0.014 to 0.020 from the 2MASS survey jHuchra et al. Il2005h . 
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Fig. 15. — CTA 1. Greyscale is Psoo/ > with a peak brightness of 45 mJy/800" beam. Con- 
tours are broadband X- rays from ROSAT, con volved to 800" . The hatched circle represents 
the cloud described by iPineault et al. I (119971 ). and the arrow indicates the new polarized 
patch not visible in the DRAO higher resolution image. 



